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Abstract. We introduce two algebraic completely integrable analogues of the 
Mumford systems which we call hyperelliptic Prym systems, because every 
hyperelliptic Prym variety appears as a fiber of their momentum map. As 
an application we show that the general fiber of the momentum map of the 
periodic Volterra lattice 

di = ai{ai_i — a^+i), i = l,...,n, a„+i = ai, 

is an affine part of a hyperelliptic Prym variety, obtained by removing n trans- 
lates of the theta divisor, and we conclude that this integrable system is alge- 
braic completely integrable. 
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1. Introduction 

In this paper we introduce two algebraic completely integrable (a.c.i.) systems, 



similar to the even and odd Mumford systems (see |12 for the odd systems and 



15 1 for the even systems). By a.c.i. we mean that the general level set of the 
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momentum map is isomorphic to an affine part of an Abelian variety and that the 
integrable flows are Unearized by this isomorphism ([0). The phase space of these 
systems is described by triplets of polynomials {u{x),v{x),w{x)), as in the case of 
the Mumford system, but now we have the extra constraints that m, w are even and 
V odd for the first system (the "odd" case), and with the opposite parities for the 
other system (the "even" case). We show that in the odd case the general fiber of 
the momentum map is an affine part of a Prym variety, obtained by removing three 
translates of its theta divisor, while in the even case the general fiber has two affine 
parts of the above form. We call these system the odd and the even hyperelliptic 
Prym system because every hyperelliptic Prym variety (more precisely an affine 
part of it) appears as the fiber of their momentum map. Thus we find the same 
universality as in the Mumford system: in the latter every hyperelliptic Jacobian 
appears as the fiber of its momentum map. 

To show that the hyperelliptic Prym systems are a.c.i. we exhibit a family of 
compatible (linear) Poisson structures, making these systems multi-Hamiltonian. 
These structures are not just restrictions of the Poisson structures on the Mumford 
system. Rather they can be identified as follows: the hyperelliptic Prym systems 
are fixed point varieties of a Poisson involution (with respect to certain Poisson 
structures of the Mumford system) and we prove a general proposition stating that 
such a subvariety always inherits a Poisson structure (Proposition ^ ). 

As an application the algebraic geometry and the Hamiltonian structure of the 
periodic Volterra lattice 

(1) fli = ai(ai_i - flj+i) i = a„+i = ai. 

Although systems of this form go back to Volterra's work on population dynamics 
( |20[] ), they first appear (in an equivalent form) in the modern theory of integrable 
system in the pioneer work of Kac and van Moerbeke (fl^), who constructed this 
system as a discretization of the Korteweg-de Vries equation and who discovered 
its integrability. Though those authors only considered the non-periodic case, we 
shall refer to (|l|) as the n-body KM system. In the second part of the paper we give 
a precise description of the fibers of the momentum map of the KM systems and 
we prove their algebraic complete integrability. 
We can summarize our results as follows: 

Theorem . Denote be M', V, T, /C the phase spaces of the (even) Mumford sys- 
tem, the hyperelliptic Prym system (odd or even), the (periodic) sl Toda lattice 
and the (periodic) KM system. Then there exists a commutative diagram of a.c.i. 
systems 



where the horizontal maps are morphisms of integrable systems, and the vertical 
maps correspond to a Dirac type reduction. 

We stress that the vertical arrows are natural inclusion maps exhibiting for both 
spaces the subspace as fixed points varieties, but they are not Poisson maps. On 
the other hand, the horizontal arrows are injective maps that map every fiber of the 
momentum map on the left injectively into (but not onto) a fiber of the momentum 
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map on the right. In order to make these into morphisms of integrable systems, we 
construct a pencil of quadratic brackets making Toda — > Mumford a Poisson map. 
For one bracket in this pencil the induced map for the KM systems is also Poisson, 
so it follows that the diagram is also valid in the Poisson category. 

A description of the general fiber of the momentum map of the KM systems as 
an affine part of a hyperelliptic Prym variety follows. Since the flows of the KM 
systems are restrictions of certain linear flows of the Toda lattices this enables us 
to show that the KM systems are a.c.i.; moreover the map $ leads to an explicit 
linearization of the KM systems. 

In order to determine precisely which divisors are missing from the affine varieties 
that appear in the momentum map we use Painleve analysis, since it is difficult 
to read this off from the map <I>. The result is that n (= the number of KM 
particles) translates of the theta divisor are missing from these affine parts. We 
also show that each hyperelliptic Prym variety that we get is canonically isomorphic 
to the Jacobian of a related hyperelliptic Riemann surface, which can be computed 
explicitly, thereby providing an alternative, simpler description of the geometry of 
the KM systems. 

The plan of this paper is as follows. In Section 2 we recall the definition of a 
Prym variety and specialize it to the case of a hyperelliptic Riemann surface with an 
involution (different from the hyperelliptic involution) . We show that such a Prym 
variety is canonically isomorphic to a hyperelliptic Jacobian and we use this result 
to describe the affine parts that show up in Section 3, in which the hyperelliptic 
Prym systems are introduced and in which their algebraic complete integrability is 
proved. In Section 4 we establish the precise relation between the KM systems and 
the Toda lattices and we construct the injective morphism We use it to give a 
first description of the general fiber of the momentum map of the KM systems and 
we derive its algebraic complete integrability. A more precise description of these 
fibers is given in Section 5 by using Painleve analysis. We finish the paper with 
a worked out example {n = 5) in which we find a configuration of five genus two 
curves on an Abelian surface that looks very familiar (Figure 2). 

As a final note we remark that the (periodic) KM systems have received much less 
attention than the (periodic) Toda lattices, another family of discretizations of the 
Korteweg-de Vries equation, which besides admitting a Lie algebraic generalization, 
is also interesting from the point of view of representation theory. It is only recently 
that the interest in the KM systems has revived (see e.g. jl^, and the references 
therein). We hope that the present work clarifies the connections between these 
systems and the master systems (Mumford and Prym systems). It was pointed out 
to us by Vadim Kuznetsov, that an embedding of the KM systems in the Heisenberg 
magnet was constructed by Volkov in . 

2. Hyperelliptic Prym varieties 

In this section we recall the definition of a Prym variety and specialize it to 
the case of a hyperelliptic Riemann surface F, equipped with an involution a. We 
construct an explicit isomorphism between the Prym variety of (F, a) and the Ja- 
cobian of a related hyperelliptic Riemann surface. We use this isomorphism to give 
a precise description of the affine part of the Prym variety that will appear as the 
fiber of the momentum map of an integrable system related to KM system. 
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2.1. The Prym variety of a hyperelliptic Riemann surface. Let F be a 

compact Riemann surface of genus G, equipped with an involution a with p fixed 
points. The quotient surface To- = F/cr has genus g' , with G = 2g' + p/2 — 1, 
and the quotient map F — > To- is a double covering map which is ramified at the p 
fixed points of a. We assume that g' > 0, i.e., <j is not the hyperelliptic involution 
on a hyperelliptic Riemann surface F. The group of divisors of degree on F, 
denoted by Div°(F), carries a natural equivalence relation, which is compatible 
with the group structure and which is defined by I? iff 2? is the divisor of 
zeros and poles of a meromorphic function on F. The quotient group Div*'(F)/ ~ 
is a compact complex algebraic torus (Abelian variety) of dimension G, called the 
Jacobian of F and denoted by Jac(F) (|^, Ch. 2.7), its elements are denoted by [V], 
where V £ Div°(F) and we write (8) for the group operation in Jac(F). Notice that 
(T induces an involution on Div°(F) and hence on Jac(F); we use the same notation 
(7 for these involutions. 

Definition 2.1. The Prym variety of (F, cr) is the (G — g')-dimensional subtorus 
of Jac(F) given by 

Prym(F/F^) ^ {[D - a{V)] \ V G Div"(F)}. 

We will be interested in the case in which F is the Riemann surface of a hyper- 
elliptic curve F^''^ : — f{x), where / is a monic even polynomial of degree 2n 
without multiple roots (in particular is not a root of /), so that the curve is non- 
singular. The Riemann surface F has genus G = n — 1 and is obtained from F'"' by 
adding two points, which are denoted by c»i and 002. The two points of F^^^ for 
which X = are denoted by Oi and O2 ■ The 2n Weierstrass points of F (the points 
(x,y) of F*^"^ for which y — 0) come in pairs (X, 0) and (— X, 0); fixing some order 
we denote them by Wi — {Xi,Q) and —Wi — (— Xi,0), where i — 1, . . . ,n. The 
Riemann surface F admits a group of order four of involutions, whose action on 
F'^^) and on the Weierstrass points {Xi,Qi) and whose fixed point set are described 
in Table 1 for n odd, n — 2g + 1 and in Table 2 for n even, n 2^ + 2 (j is the 
hyperelliptic involution) . 



Table 1: n odd 







Oi 


O2 




002 


w^ 


-Wi 


Fix 


I 


{x, -y) 


O2 


Oi 


002 


OOl 


w. 


-w. 


W„ -Wi 


a 


{~x,y) 


Oi 


02 


002 


OOl 


-w. 


w. 


Oi, O2 


T 


{-X, -y) 


O2 






002 


~w. 


w. 


OOl, 002 



Table 2: n even 





(x,y) 


Oi 


O2 


OOl 


002 


w. 


-w. 


Fix 


I 


{x, -y) 


O2 


Oi 


002 


OOl 


w. 


-w. 


w„ -w. 


a 


{-x,-y) 


O2 


Oi 


002 


OOl 


-w. 


w. 




T 


i-x,y) 


Oi 


02 


OOl 


002 


-w^ 


W^ 


Oi, O2, OOl, 002 
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For future use we also point out that for points P G F which are not indicated on 
these tables, neither cr(P) nor r(P) coincide with i{P). 

The involutions a and r lead to two quotient Riemann surfaces F^^ := F/cr and 
F^ := F/r, and to two covering maps n„ : T and tt,- : F ^ F^. It follows from 

Tables 1 and 2 that the genus of F,- equals g, while the genus g' of F^- is g or g + 1 
depending on whether n is odd or even. Also, the dimension of Prym(F/Fcr) = g 
(whether n is odd or even). If the equation of F^^^ is written as y"^ = g{x'^) then 
for n odd, F^"' has an equation = g{u) while F^"-* has an equation ~ ug(u); 

for n even the roles of F^"'' and Fr"' are interchanged. 

In order to describe Prym(F/Fcr), which we will call a hyperelliptic Prym variety, 
we need the following classical results about hyperelliptic Riemann surfaces and 
their Jacobians (for proofs, see Ch. Ilia). 

Lemma 2.2. Let V be a divisor of degree H > G onT, where G is the genus ofV, 
and let P he any point on F. There exists an effective divisor £ of degree G on T 
such that 

V^£ + {H- G)P. 
Corollary 2.3. For any fixed divisor Vq of degree G, Jac(F) is given by 
Jac(F) = j ^P-Po l^^eFl. 



Lemma 2.4. Let V be a divisor on F of the form T) = X^iLi {Pi ~ Qj) where H < G 
and Pi ^ Qj for all i and j . Then [D\ — if and only if H is even and V is of the 
form 

H/2 

i=l 

for some points Ri, Si G F. 

2.2. Hyperelliptic Prym varieties as Jacobians. In the following theorem we 
show that for any n the Prym variety Prym(F/Fo.) associated with the hyperelliptic 
Riemann surface F is canonically isomorphic to the Jacobian of Fr. 

This result was first proven by D. Mumford (see ||l^) for the case in which 
TTcr : F ^ Fcr is unramificd (n even) and by S. Dalaljan (see 0) for the case in 
which TT^ : F ^ F^ has two ramification points (n odd). Our proof, which is valid 
in both cases, is different and has the advantage of allowing us to describe explicitly 
the affine parts of the Prym varieties that we will encounter as affine parts of the 
corresponding Jacobians. 

Theorem 2.5. Let tt* denote the homomorphism Div'^(FT-) — > Div'^(F) which sends 
every point ofT^ to the divisor on F which consists of its two antecedents (under t). 
The induced map 



n : Jac(F^) ^ Prym(F/F„) 
[V] ^ [<P] 



is an isomorphism. 
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Proof. It is clear that the homomorphism 11 is a weh-defined: if [D] — then V 
is the divisor of zeros and poles of a meromorphic function / on , hence tt*!? is 
the divisor of zeros and poles of / o r and [n*V] = 0. To see that the image of 11 
is contained in Prym(r/rcr), just notice that tt*{'D) can be written as £ + t{£) for 
some £ € Div"(r), so that 

[<(P)] ^[£ + t{£)] = [£- a{£)] G Prym(r/r.). 

Since Jac(rT-) and Prym(r/rcr) both have dimension g it suffices to show that 11 is 
injective. Suppose that [vr*!?] = for some T ) £ Div"(rT-). We need to show that 
this implies [D] = 0. It follows from Corollary 2.3 that we may assume that V is of 

2gcx)i 
= ooi, 

□ 



the {orinJ2i=iPi-9T^T{coi), where Pi G T^. Then Tr*D = yj _-, Pj+TjP.j) 



{Tr^{Pi) = Pi). Since 2g < G and ooi ^ «(ooi) Lemma 2.4 implies that Pi 
i.e., Pi — 7r^(ooi) for all i. 

2.3. The theta divisor. We introduce two divisors on Jac(r) by 



(2) 



(3) 



62 = 



0O2 



Gooi 



P^eT 



P^eT 



These two divisors are both translates of the theta divisor and they differ by a shift 
over [002 — ooi]. Since 002 = «(c!Oi) they are tangent along their intersection locus, 
which is given by 



n = 



■G-2 

E 



Pi + 002 - (G - l)cx)i 



P^ er 



Proposition 2.6. When n is odd Prym(r/rcr) fl (Oi U02) consists of three trans- 
lates of the theta divisor 0/ Jac(r^), intersecting as in the following figure. 




[001 + O] 



[002 + O] 



Figure 1 

Proof. We use the isomorphism 11 to determine which divisors of Jac(rT) get 
mapped into Oi and 02- Since Oi and O2 are the only points of F on which i 



and T coincide. Lemma 2.4 implies that the only divisors D — ^f^iPi — giTriooi) G 



Div(r^) for which 7r*P contains, up to linear equivalence, cx^i or 002 are those for 
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which at least one of them contains 7rT-(ooi) or 7rT-(oo2) or nr{Oi) {^'^t{02))- De- 
noting O — TTriOi) we find that these points constitute the foUowing three divisors 
on Jac(r^). 

Oi^l - (g - l)7r^(ooi) 

I Li=i 

9^1 y^pi + O - ff7rr(0Ol) 




They all pass through 



3-2 



^Pi + 7rr(oo2) - (5 - l)7r^(ooi) 



which is the tangency locus of 9i and 92, and 9i intersects 9 in addition in 

LOi^l ^p.i + TTrioOi) + O - giTrioOi) | G L 

which is a translate of w. □ 
When n is even then clearly Prym(r/ro-) is contained in Oi, but the following 
result, similar to Proposition holds for an appropriate translate of Prym(r/rcr). 
The proof is left to the reader. 

Proposition 2.7. When n is even and i G {1,2} then (Prym(r/rCT) ^ [Oi — ooi]) n 
(01 U 62) consists of three translates of the theta divisor 0/ Jac(rT-), intersecting 
as in Figure 1 (in which O should now be replaced by 02)- 

We will see in the next section how in both cases (n odd/even) the afRne 
variety obtained by removing these three translates from the theta divisor from 
Prym(r/r£r) can be described by simple, explicit equations. 

3. The hyperelliptic Prym systems 

In this section we introduce two families of integrable systems, whose mem- 
bers we call the odd and the even hyperelliptic Prym systems, where the adjective 
"odd/even" refers to the parity of n, as in the previous section, and where "hy- 
perelliptic Prym" refers to the fact that the fibers of the momentum map of these 
systems are precisely the affine parts of the hyperelliptic Prym varieties that were 
considered in the previous section. These systems are intimately related to the even 
Mumford systems, constructed by the second author (see ffaj), as even analogs of 
the (odd) Mumford systems, constructed by Mumford (see [Q). 

3.1. The Mumford systems. We first recall the definition of the g-dimensional 
odd and even Mumford systems and we describe their geometry. Details, general- 
izations and applications can be found in |16| . The phase space of each of these 
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systems is an affine space C , which is most naturally described as an afRne space 
of triples {u{x),v{x),w{x)) of polynomials, often represented as Lax operators 

\ u{x) —v{x^ 

where u{x), v{x) and are subject to certain constraints. Denoting by Mg 

(resp. -M'g) the phase space of the g-th odd (resp. even) Mumford system these 
constraints are indicated in the following table. 



Table 3 



phase space 


dim 


u{x) 


v{x) 


w{x) 


Mg 




monic 




monic 


35 + 1 




deg < g 








deg = g 




deg = g + 1 


M'g 


3g + 2 


monic 
deg = g 


deg < g 


monic 
deg = c/ + 2 



It is natural to use the coefficients of the three polynomials u{x),v{x),'w{x) as 
coordinates on Mg and on M'g: for M'g for example, which will be most important 
for this paper, we write 

u{x) =x^ + Ug-ix^~^ H h Wo, 

v{x) = Vg-ix^~'^ H \-vo, 

w{x) = a;S+2 + Wg+ix^+'^ + --- + wo, 

or, in terms of the Lax operator L{x), as 




E 



Vi Wi 



X . 



/ 0<i<g \ 



We define on Mg and on M'g a map H with values in (a finite-dimensional aflane 
subspace of) C[a;] by 

H{L{x)) = - dct{L(x)) = u{x)w{x) + v^{x). 

Notice that the degree of H, as a polynomial in x, is odd for Mg and even for 
M'g, explaining the terminology odd/even Mumford system. Writing H{L{x)) = 
2;2s+i _j_ H^gX^^ + ■ • • + Hq in the odd case, and similarly in the even case, wc can 
view the Hi as functions on phase space; clearly they are polynomials in the phase 
variables Ui, Vi, Wi. These functions are independent and they arc in involution with 
respect to a whole family of compatible Poisson brackets on Mg (or M'g). Namely, 
it is easy to check by direct computation that for any polynomial (f of degree at 
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most g the foUowing formulas 

{u{x),u{x')}'^ = {v{x),v{x')}'^ = 0, 

(4) {uix),w{x')} 



^ x-x' 



X ~ X 

{v{x), w(x')Yj^ = '^^ / - a{x + x')uixMx'), 

X ~ x' 

{w{x),w{x')Yj^ = 2a{x + x') {v{x)(p{x') ~ v{x')(p{x)) , 

define a Poisson structure on A4g and on M'g, where a — 1 for the odd case (i.e., 
on Aig) and a ^ x + Wg+i — for the even case (i.e., on A4'g)', we use the same 
notation {•, -j^ for the Poisson brackets on both spaces. When (fi then the 
rank of this Poisson structure is 2g on a Zariski open subset. Since 

(5) {H{x), H{x')Ym = {u{x)w{x) + v\x), u{x')w{x') + v\x')}l^ = 0, 

the functions Hi are in invohition. For any <y9 ^ of degree at most g and for any 
y (z C which is not a root of ip, let 

Then 

1 

is a Hamiltonian vector field, independent of the choice of the Poisson structure 
{•, - j^, i.e., it is a multi-Hamiltonian vector field. Explicitly, Xy can be written as 
the following Lax equation 

(7) XyL{x) = -^[L{x),L{y) + (x - y)B{x, y% 

X y 

where 

a{x + y)u{y) 



(6) Xy = {-, H^{y)VM = -7-t{- > ^ivMy) + v'iy)VM 



B(x,y) = 







In view of (^) the vector fields Xy corresponding to any two values of y commute 
and a dimension count shows that the Mumford systems are integrable in the sense 
of Liouville. 

In fact, the Mumford systems are algebraic completely integrable systems (a.c.i. 
systems), meaning that the general level of the momentum map H (equivalently, 
W^) is an affine part of an Abelian variety and that the fiow of the integrable vector 
fields Xy is linear on each of these Abelian varieties. A precise description of these 
fibers is given in the following proposition. 

Proposition 3.1. Let f{x) be a monic polynomial of degree 2g + 1 (resp. 2g + 2) 
without multiple roots and let T denote the Riemann surface corresponding to the 
(smooth) affine curve (of genus g) defined by F*-*^) : y^ — f{x). Then the fiber over 
f{x) of H : A4g — > C[x] (resp. H : M'^ — > is isomorphic to Jac(r) minus 

its theta divisor (resp. minus two translates of its theta divisor which are tangent 
along their intersection locus). 
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Proof. We shortly indicate the idea of the proof, since it wiU be useful later. To 
v{x), w{x)) in the fiber over /(x) of : A^g C[a;] one associates a divisor 
V — Yli=ii^i^yi) ~ r (where {00} := F \ F^"^) by taking for Xi the roots 

of u{x) and yi = v{xi). This map is injective and its image consists of those 
divisors ~ 5°° ^'^^ which Pi e F^"^ and for which i j ^ Pi ^ ^Pj)- 

Mapping V to its equivalence class [D] we get an injective map into Jac(F) and 



the complement of its image is the theta divisor | X]i=i Pi ^ (.9 ^ l)oo 
For {u{x),v{x),w{x)) G A4'g the construction is similar but the complement of the 
image consists of two translates of the theta divisor because F \ F*^"^ consist now of 
two points. □ 
Notice that the fact that the general fiber of the momentum map is g-dimensional 
implies our earlier assertion that the functions Hi are independent, an essential 
ingredient in the proof of the Liouville intcgrability of the Mumford systems. 

3.2. Phase space and momentum map. In this paragraph we introduce the 
phase spaces of the hypcrelliptic Prym systems, which can be seen as affine sub- 
spaces of the even Mumford systems; they are not Poisson subspaces with respect 
to the Poisson structures (Q), as we will discuss in the next paragraph. However, the 
moment map will just be the restriction of the moment map H to these subspaces. 

For n > 1, let us denote by Vn and by V'^ the set of triples of polynomials 
{u(x),v{x),w{x)) which satsify the constraints indicated in the following table: 



Table 4 



phase space 


dim 


u(x) 


v{x) 






Pn 


3n + 1 


monic 


odd 


monic 




even 
deg = 2n 


deg < 2n 


even 

deg = 2n 4 


-2 






monic 




monic 




' n 


3n + 2 


odd 
dcg = 2n+l 


even 
deg < 2n + 1 


odd 

deg = 2n- 


-3 



Comparing this table to Table 3 one sees that Vn is an affine subspace of A^2„ 
and that Vl^ is an affine subspace of Aljn+i- For reasons that will become clear 
shortly, we will say that Vn (resp. V'^) is the phase space of the n-th odd (resp. 
even) hyperelliptic Prym system. We denote the restriction of H : A^2„ ^ C[x] 
to Vn as well as the restriction of H : M.'2n+\ ~^ 'C\x\ to V[^ also by PI. Notice 
that iJ, which is on these subspaces still given by 

H{L{x)) = -detL{x) = u{x)w{x) +v'^{x), 

takes now values in C[x^]. It follows that the corresponding Riemann surface F is 
of the type considered in Section We show in the following two propositions that 
the fibers of H are affine parts of hyperelliptic Prym varieties. 

Proposition 3.2. Let f{x) be a monic even polynomial of degree Ag + 2 without 
multiple roots and let F denote the Riemann surface corresponding to the (smooth) 
affine curve (of genus 2g) defined by T^^^ : y^ — f{x). The fiber of H :Vg ^ C[a;^] 
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over f{x) is isomorphic to Pryni(r/rcr) = Jac(rT-) minus three translates of its 
theta divisor, intersecting as in Figure 1. 

Proof. Since the fiber over f{x) of H : Vg ^ C[x^] is contained in the fiber over f{x) 
of H : M'2g C[x] it is a subset of Jac(r). In fact it is a subset of Prym(r/rCT). 

To see this, consider the divisor D = X]?=i (■^i' J/O ^ Syooi which corresponds to a 
triple (u{x),v{x),w{x)), with u,w even and v odd. The roots of u come in pairs 
{xi,Xj = —Xi) and yj = v{xj) = v{—Xi) = —v{xi) = —tji (recall that can never 
be a root of u because then / would have as a double root), hence the points in T) 
come in pairs P,t{P) and [D] belongs to Prym(r/rCT). The points of Prym(r/rcr) 

which do not belong to the fiber are those E^=i^«-2gooi e Prym(r/r„) for 

which at least one of the Pi equals ooi or 002, i.e., the points on Oi U 02, as 
defined in (||) and By Proposition 2J the fiber is isomorphic to an affine part 



of Prym(r/rcr) obtained by removing three translates of the theta divisor. □ 

Proposition 3.3. Let f{x) be a monic even polynomial of degree 4g + 4 without 
multiple roots and let T denote the Riemann surface corresponding to the (smooth) 
affine curve (of genus 2g A- 1) defined by V'^^^ : y^ ~ f{x). The fiber over f{x) 
of H : V'g ^ C[x^] is reducible and each of its two components is isomorphic to 
Prym(r/r(j) = Jac(rr) minus three translates of its theta divisor, intersecting as 
in Figure 1. 

Proof. Consider the divisor T) — Y^i^^ i^iTVi) ~ (^5 + l)oo which corresponds to 
a triple {u{x), v{x), w{x)), with m, w odd and v even. is a root of u and its other 
roots come in pairs Xi,Xj = —Xi and yj = v{xj) = v{xi) = yi, hence the points in 
V consist of Oi or O2 and the others come in pairs P,t{P). It follows that [D] 
belongs to Prym(r/ro-) ® [Oi - ooi] or to Prym(r/rCT) ® [Oi - 002]. The points of 

Prym(r/rcr) which do not belong to the fiber are those Pi " (^5 + l)ooi 

for which at least one of the Pi equals ooi or 002, i.e., the points on Qi U 02, 
as defined in (||) and (||). By Proposition 2_i the fiber consists of two copies of 



an affine part of Prym(r/rcr) obtained by removing three translates of the theta 
divisor. Notice that the fact that the fiber is reducible can also be deduced from 
the fact that /(O) = u(0)w(0) + ^^(O) = wg. □ 

3.3. Flovifs and Hamiltonian structure. We now show how the brackets (jj) 
lead to a family of compatible Poisson structures on Vn and on V^^, we construct 
the integrable vector fields for the hyperelliptic Prym systems and we establish 
their algebraic complete integrability. Consider the natural inclusions Pg '-^ ■^2g 
and P'g ^ A^2g+i consider the involution j : Ai',^ Ai'^ defined by: 



v{x) w{x) \ I —v{—x) w{—x) 
u{x) —v{x) I \ u(~x) v{—x) 



Then we see that the image of Pg ^ -^2g the fixed point variety of j, while the 
image of V'g ^ .^2^+1 is the fixed point variety of —j. We claim that j (resp. —j) 
is a Poisson automorphism of {M2g, {■, ■}%) (resp. (A^2g+i, {•, whenever ip 
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is an even (resp. odd) polynomial. In fact, taking ip even we have 

{u(x),v{x')Ym 



{u{x) o J, v{x') o j}^ = {u{~x), ~v{~x')}'^ 

u{—x)(f{x') — u{~x')ip{x) 



X — x 

showing that for any i and j, {ui o j,Vj ° j}^ = {ui,Vj}'^ and similarly for the 
Poisson brackets of the other components. Since all brackets are linear in ip the 
result for ip odd also follows, when j is replaced by —j. 

The following proposition, which can be seen as a particular case of Dirac reduc- 
tion, yields a Poisson structure on the fixed point variety of a Poisson involution. 
For the general theorem on Dirac reduction, see Weinstein (|^, Prop. 1.4) and 
Courant (Q, Thm. 3.2.1). For our convenience we give a proof in the algebraic 
category; the proof is easily adapted to smooth manifolds. 

Proposition 3.4. Suppose that (M, {•,•}) is an affine Poisson variety, equipped 
with an involution j which is a Poisson map. Let N be the subvariety of M con- 
sisting of the fixed points of j and denote the inclusion map N ^ M by i. Then N 
carries a (unique) Poisson structure {•, -j^ such that 

(8) i*{F,G}^{i*F,i*G}^ 

for all F,G ^ 0{M) that are j-invariant. 

Proof For f,g G 0{N) we choose F, G G 0{M) such that / = i*F and g = i*G. 
We may assume that F and G are j-invariant by replacing F by {F + j*{F))/2 
and similarly for G. We define {/, g}j^ = i* {F^ G} and show that this definition is 
independent of the choice of F and G. To do this it is sufficient to show that if G is 
J-invariant and i*F — 0, then i* {F, G} — 0. Since the ideal of functions vanishing 
on N is generated by j-anti-invariant functions {]*F = —F) it suffices to show this 
for F G 0{M) such that fF = —F. By assumption j is a Poisson map, ] o i — ] 
and f'G — Gso that 

I* {F, G] = i*f {F, G} = I* {fF, fG} = -I* {F, G} , 

showing our claim. Similarly, the bracket of any two j-invariant functions is j- 
invariant. In view of this and because the definition of {•, - j^ is independent of the 
choice of F and G we have for any f,g,h G 0{N) that 

{{f,g}j„h}^^t*{{F,G},H}, 

leading at once to the Jacobi identity for {•, Similarly the fact that {•, - j^ is 
an anti-symmetric biderivation follows. □ 
The Hamiltonian structure of the hyperelliptic Prym systems and its algebraic 
complete integrability is described in the following proposition. 

Proposition 3.5. Let if be an even (resp. odd) polynomial of degree at most 2g + l, 
(p ^ 0. The Poisson bracket {•, - j^ on M.'2g, (resp. on ■M'2g+i)j given by induces 
a Poisson bracket j-,-}^ on Vg (resp. on Vg), the components of H :Vg ^ 'C[x'^] 
(resp. H : V'g ^ ^[x'^]) are in involution and they define an (algebraic completely) 
integrable system on Vg (resp. on Vg). 

Proof. Since for ip even (resp. odd) the image of Vg ^ -^2g (resp. Vg ^ -^2g) 
is th e fix ed point set of the Poisson involution ] (resp. -]) it follows from Proposi- 
tion 3^ that Vg (resp. V'g) inherits a Poisson bracket from A^2g, which we denote 
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in both cases by {•, - jp. We exemplify the computation of the reduced brackets by 
deriving the formula for {u{x), on Pg {ip even). Notice that since u is even 

and V is odd the polynomial {u{x),v{x)}'^, which is a generating function for the 
Poisson brackets {ui^Vj}"^, is even in x and odd in x' . Obvious j-invariant exten- 
sions of the functions U2i and U2i+i are the corresponding functions U2i and U2i+i 
on ■M'2g- Therefore {u{x),v{x')}p is computed by taking in {u{x),v{x')}'^ the 
terms that are even in x and odd in x' and restricting the resulting polynomial to 
the image of Vg, as embedded in A^2g- Using the fact that the terms of a bivariate 
polynomial F(x,x') that are even in x and odd in x' are picked by taking 



we find for 



- x') + F{~x, x') - F{x 



F{x,x') = {u{x),v{x')}'^ 



')~Fi-x,-x')) 
u{x)ip{x') — u{x')(p{x) 



that the reduced Poisson bracket, for ip even, is given by 



{u{x),v{x')}^ = — 



x' f {u{x) + u{—x))ip{x') — {u{x') + u{~x'))ip{x) 



,u{x)(p{x') — u{x')ip{x) 



x^ — x'^ 



Repeating the same computation for the other coordinates we find the following 
formulas for {•, -jp, 

= {v{x),v(x')Y^ = 0, 

iu(x)Lp(x') — u(a;')(p(a;) 



{u{x),v{x')) 
{u{x),w{x')\ 
{v{x\ w[x'')) 



= -2 



a;^ — X 
xv{x)(p{x') — x'v{x')(p{x) 



,.'2 



X^ — X 

w{x)ip{x') — w{x')ip(x) 



xu{x)(p{x'), 



x-^ — X' 

{w{x), — 2 {xv{x)lp{x') — x'v{x')(p{x)) . 

Using the fact that all Poisson brackets are linear in (p one finds that the formulas 
for the reduced bracket on V'g (with (p odd) are formally identical to the above 
ones. It is now obvious that the components of the new momentum map H are 
in involution. Since we know that the fibers of H : Vg ^ C[x'^] are afhne parts of 
Abelian varieties of dimension g, the components of the new H are independent. 
The integrable vector fields Xy on Vg are computed from {• , H{y)}^, to wit 



(9) XyL{x) = 



1 



y 



L{x), 



yv{y) xw{y) + x{x'^ - y'^)u{y) 

xu{y) -yv{y) 



Since the formulas for the reduced brackets on Vn and on V'g are formally the same 
the vector fields Xy on V'g are also given by @. Finally, the flows of all vector 
fields Xy are linear since they are restrictions of linear flows, showing that the odd 
hyperelliptic Prym systems are algebraic completely integrable. □ 
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4. The periodic Toda lattices and KM systems 



In this section we show that the (periodic) KM systems are related to the (peri- 
odic) si Toda lattices in the same way as the hypcrclliptic Prym systems arc related 
to the even Mumford systems and we construct a morphism from the Toda lattices 
to the even Mumford systems, which induces a morphism from the KM systems to 
the odd or the even hypcrclliptic Prym systems. The latter map is then used to 
describe the level sets of the momentum map of the KM systems. 

4.1. Prom Toda to KM. The phase space of the periodic sl(n) Toda lattice (n- 
body Toda lattice for short) is the afRne variety of all Lax operators in si{n)[h, h~^] 
of the form 

/ 



(10) 



L{h) 



bi 
1 



ai 











62 a2 







V 



bn-l Cln-l 

httn 1 bn J 

with niLi oij = 1- It carries a natural Z/n action, defined by 

((ai, a2 . . . , an), (61, 62 • • • , &«)) ((012, as • • • , oi), (62. bs, . . . , 61)). 
For this reason it is convenient to view the indices as elements of Z/n and we put 



do, bn 



bo, fln+i = ai, 6„+i = 61, . . . Define 
1 



1 + i 



,tr{L{hr+'), 



0, ...,n-2. 



In- 



-tr(L(/i)")-/i-i 
n n 



and notice that these n functions are independent of h. They are in involution with 
respect to the linear Poisson bracket {•, -j^, defined by {ai,aj}^ = {bi,bj}}^ = 0, 
{ui, bj}]j- = ai{5ij — 5i+ij), which has /q as a Casimir. They are also in involution 
with respect to the quadratic Poisson bracket {•, -j^, defined by 



{ai,6j}^ = aibj{5^j - (5i+ij), 



{bi,bj}^ = ai{Sij+i — Si+ij), 

which has deti as a Casimir. Since j-,-}^ and {•, -j^ are compatible we may 
define, for any ip € iC[x] of degree at most 1 a Poisson bracket on 7^ by {•, •} 



r 



{■> "It + fo {■, -It) where (p{x) = (pix + (po- 
The commuting vector fields Xi = {■ , lijlj- admit the Lax representation 

(11) XMh) = [m,{Lih)%], 

where the subscript + denotes projection into the Lie subalgebra of sl{n)[h,h~^] 
generated by the positive roots, i.e., 

(12) =E^^^'+^u(^o)' 

i>0 
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where su(j4o) denotes the strictly upper triangular part of Aq. The vector fields 



Xi are also Hamiltonian with respect to {• 
general fiber of the momentum map K : Tn 



and their flows are linear on the 
-> C[x], which is defined by 

det{xld -L{h)) = -h- ^ + K{x)/2; 

since the general fiber of K ia an afhne part of a hyperelliptic Jacobian, the n-body 
Toda lattice is an a.c.i. system (see Q for details). For higher order brackets for 
the Toda lattices, see 

We now turn to the n-body, periodic, Kac-van Moerbeke system (n-body KM 
system, for short). Its phase space is the subspace of %i consisting of all Lax 
operators ( |lo|) with zeros on the diagonal. /C„ is not a Poisson subspace of 7^. 
However, /C„ is the fixed manifold of the involution j : Tn ^ Tn defined by 

((oi, a2 . . . , a„), (61, 62 . . . , 6„)) ((ai, 02 . . . , a„), {~bi, -62 ■ • • , -&«)), 



which is a Poisson automorphism of (7^, {•, -j!^). Therefore, by Theorem 3.4, /C„ 
inherits a Poisson bracket {•, - j^ from {•, ■}^, which is given by 

{ai,aj}j^ = aiaj{dij-^i — Si+ij). 

It follows that the restriction of the momentum map K to /C„ is a momentum map 
for the n-body KM system. Notice that Ij = for even j, while for j odd the Lax 
equations ( pi] ) lead to Lax equations for the n-body KM system, merely by putting 
bi = ■ ■ ■ = bn = 0. Taking j = 1 wc find the vector field 

(13) fli = ai(aj_i - a.j+i), i = l,...,n, 

which was already mentioned in the introduction. More generally, taking j odd we 
find a family of commuting Hamiltonian vector fields on /C„ which are restrictions 
of the Toda vector fields, while for j even the Toda vector fields Xj are not tangent 
to ICn- In order to conclude that the KM systems are a.c.i. we need to describe 
the fibers of the momentum map K : /C„ — » C[x]. This will be done in the next 
paragraph. 

4.2. Algebraic integrability of KM. We first define a map ^ : Tn ^ -^n-i 
which maps the n-body Toda system to the even Mumford system. The following 
identity, valid for tridiagonal matrices, will be needed. 

Lemma 4.1. Let M be a tridiagonal matrix, 

( R. rv. f\ ... n \ 



M 



\ 



Pi 


ai 


• 








71 


P2 


a2 










72 















Pn-l 












■ 7n-l 





J 



and denote by A^j^.. 
from M the rows ii, 



the determinant of the minor of M obtained by removing 
, ik and the columns ii, . . . , z^. Then: 



(14) 



AiA„ - AAi, 
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Proof. For n ~ 2 this is obvious. For n > 2 one proceeds by induction, using the 
foUowing formula for calculating the determinant A of M, 

(15) A = /3„A„ - q;„_i7„_iA„_i^„. 

□ 

In the sequel we use the notation Aij^. . ^^. from the above lemma taking as M the 
tridiagonal matrix obtained from xld —L{h) in the obvious way, i.e., by removing 
the two terms that depend on h. In this notation the characteristic polynomial of 
L(h) is given by 

(16) det(a; Id -L{h)) ^ -h - h^^ + A - a„ Ai,„. 
Proposition 4.2. For any m = 1, . . . , n the map : Tn -^'n-i defined by 

U{x) = Am, 

(17) v{x) 

w{x) = (a; - 6„i)^A„i + 2{x - 6,„)(a„i_i A„i_i_,„ + amA„i,m+i) 

maps each fiber of the momentum map K : Tn C[x] into a fiber of the momentum 
map H : AA'^-i — > C[a;]. The restriction of to /C„ takes values in V n-i when n 
is odd and in Vn^^ when n is even, mapping in both case the fiber of the momentum 
map K : ICn C[a:] into the fiber of the momentum map H : V n-i — > C[a;^] (or 
H : V'n__^ C[a;^]j. As a consequence the general fiber of the momentum map of 
the KM systems is an affine part of a hyperelliptic Jacobian. 

Proof. Since the momentum map is cquivariant with respect to the "L/n action on 
Tn it suffices to prove the proposition for m — n. 

It is easy to see that the triple (u, v,w), defined by ( p7| ) satisfies the constraints 
u, w monic, deg w = degu + 2 = n + 1 and degv < n — 1, so that $„ takes values 
in A4'j_]^. Moreover, taking /3i = ■ • ■ = /3„ = a; in ( p^ ) implies that when all entries 
on the diagonal of L{h) are zero then A^^^ has the same parity as rt — p, so that 
the triples {u,v,w) which correspond to points in /C„ have the additional property 
that V has the same parity as n while u and w have the opposite parity. Therefore 
the restriction of to ICn takes values in V n-i when n is odd and in Vn._, when 
n is even. 

For p{x) a monic polynomial of degree n, let L{h) G K^^{2p{x)), i.e., 

(18) p{x) = (x - 6„)A„ - a„Ai„ - a„_iA„„i_„. 

Proving that $„(L(/i)) belongs to iJ^^(p^(x)— 4) amounts to showing that u{x)w{x) + 
v'^{x) = p'^{x) — 4, which follows from a direct computation, using ([ij). The com- 
mutativity of the following diagram follows: 

Tn ^M'_, 



K 



H 



where (j) is defined by (j){q) — {q/2)'^ — 4, for q G C[a;]. 
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To show that the map <i>„ is injective let {u{x),v{x),w{x)) e $„(7^i). We show 
that the matrix L(h) € 7^ which is mapped to this point is unique. 

First observe that the monic polynomial p{x) = A — a„Ai.„ can be recovered 
from u{x)w{x) + v{x)^ = p{x)^ — 4. We can then determine 6„ from the following 
two formulas: 

p{x) = - 

uix) - A„ = - 6, j + . . . . 

Next, the second relation in ( p7| ) and ( p^ lead to the system: 

2„_iA„_i,„ - a„Ai,„ = 
2„_i A„_i,„ + a„Ai,„ = (x - 6„)?i(a;) - p{x). 

This linear system completely determines the products a„Ai^„ and a„_iA„_i^„. 
Because the determinants of the principal minors of xld —L{h) are monic poly- 
nomials, this means that we know a„, Ai „ and A„_i „ separately. From A — 
p{x) + a„Ai_„ we also obtain A. 

We have now shown how 5„, a„. A, A„ and A„_i,„ are determined. We proceed 
by induction, showing how to determine bn-k-i, o-n-k-i, ^n-k-i,...,n once we 
know bn~i, Un-i and A„_i^...^„ for z = 0, . . . , fc. We use (^5|) to obtain the recursive 
relation: 

^n-k+l,...,n = (a; — bn-k)^n-k,...,n ~ 0,n-k-l^n-k-l,...,n- 

This determines the product an-k-i^n-k-i,...,ni but also a„_fe_i and A„_fe_i^...^„ 
separately, again because An-k-i,...,n is monic. Now from An-k-i,...,n and An-k,...,n 
we know, as above, the sums X]r=i^ ^ ^^'^ "^7=1 ^ Hence, bn-k-i is deter- 
mined. □ 
We saw in Proposition ^.3| that the fibers of the momentum map of the even 
Prym system are reducible (two isomorphic pieces), so there remains the question 
if the same is true for the 7i-body KM system for even n. To check that this 
is so, note that the highest degree coefficient of the characteristic polynomial of 
L(h) gives, for n even, the first integral / — aia^a^ ■ ■ ■ a„_i -I- 020405 • • • o„. Since 
ai02...o,i = 1, for generic values of /, the variety defined by 

010305 • • • o„_i — constant, 020405 • ■ ■ o„ = constant, 

is reducible, and the claim follows. Note however that both 010305 • • -On-i and 
02O4O6 • • - On are first integrals themselves, so we can construct a momentum map 
using these integrals (instead of their sum and product) and then the general fiber 
is irreducible. 

The map $m '■ %i — * A4^_i not only maps fibers to fibers of the momentum 
maps, but it maps the whole hierarchy of Toda flows to the Mumford flows defined 
by (^. To see this, we construct a family of quadratic Poisson brackets {•, -j^ ^ 
on which make this map Poisson. 

First observe that there exist unique polynomials p{x) and r(x), with p{x) monic 
of degree n and r(x) of degree less than n, such that 

(19) u{x)w{x) + v{xY — p{xY + r{x). 
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The coefiicients of p{x) and r{x) are regular functions of Ui, Vi and w;. Hence, 
we can define a skew-symmetric biderivation on the space of regular functions of 
M'n-i by setting, for any (p € C[x] of degree at most 1, 

{uix),uix')rM^^^{vix),vix')rM^O, 

{u{x),v{x')}'^ ,j = {u{x), v{x')Yj^ + a'^{x + x')u{x)u{x'), 

{u{x), = {u{x), w{x')Y^ ~ 2a'^{x + x')u{x)v{x'), 

{vix), wix')}l, ,^ = {v{x),w{x')}Z + a^{x + x')u{x)w{x')), 

{w{x),w{x')rM^^ = {w{x),w{x')}Z + 2a^{x + x') {w{x)v{x') - w{x')v{x))), 

where a^{x) = ip{a{2x)/2). Notice that the polynomial pip^ used in the definition 
of the bracket, depends on the phase variables. 

Proposition 4.3. Let ip he a polynomial of degree at most 1. Then 

(i) g Poisson bracket on A4'^^_i and the maps 

are Poisson and map the Toda flows to the Mumford flows; 

(ii) For If odd, the bracket {•, - j^ ^ induces a Poisson bracket {•, ■}-p ^ on P(n-i)/2 
(resp. on ^'^^2-1)' '^'^'^ maps 

■■ (/C„,{-,-}^) (V(n^i)/2,{-r}-p.q) 

are Poisson and map the flows of the n-body KM system to the flows of the 
hyperelliptic Prym systems. 

Proof. We take the bracket of both sides of ( |l9|) with u{x) to obtain 

My my) ^ = My) Mx),p{y)}x^^g + {u{x),r{y)}i^^g. 

X y 

ic. ^;-,r;c.;Kl<a K„ QJ^r.^ 



It follows that {u{x), r{y)}'^ ^ is divisible hy p{y). Since {u{x), r(y)}^ ^ is of degree 



less than niny and smcep{y) is monic of degree n we must have {u(x), r{y)}'^ n ~ ^ 



and 



r / N / u{x)viy) - u{y)v{x) 

{u(x),p{y)}M^g = — — p{y). 

X y 



Similarly, we find {w(x),r(y)}^ = r(j/)}^ == and also that: 



{v{x),p{y)}j^^g = ^ w{x)u{y)^ ^u{x)w{y) ^ ^ y)u{x)u{y)^ , 

r / N / M / N / v(x)w(y) — w(x)vly) , , ^ , \ 

{w{x)My)}M,q - viy) [ + "(^ + yX^Xy)\ ■ 

These expressions also allow one to compute the brackets of u(x), v{x), wlx) and 
p{x) with a{y), and the check of the Jacobi identity follows easily from it. Therefore, 
{•, - j^ g is a Poisson bracket for which the coefficients of r{x) are Casimirs. 

If we compare the expressions above for the brackets with p{y) with expressions 
for the Mumford vector fields, we conclude that they are Hamiltonian with 
respect to {■, -j^ ^ with Hamiltonian function K. Checking that is Poisson can 
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be done by a straightforward (but rather long) computation using the foUowing 
expressions for the derivatives of Aij^...^i^_: 

9Ajj^_^ _ I -Ai,i+i_i^_...,i^, + l ^ {ii, . . . ,ik} , 
1 otherwise, 

9Aij^_^ _ I -Ai,ii,...,i^, i ^ {ii, . . . ,ik} , 

I otherwise. 

For the second statement, one easily checks that when is odd then j is a Poisson 
involution, so that there is an induced bracket on 'P(„_i)/2 or on 'P^/2_i- Explicit 
formulas for this bracket are computed as in the proof of Proposition |3.5| . The other 
statements in (ii) then follow from (i). □ 

It is easy to check that the Poisson brackets {•, - j^ ^ and {•, - j^ on M'n_i are 
compatible, when and tp have degree at most 1. This is however not true when 
Tp is of higher degree. 

5. Painleve analysis 

The results in the previous section show that the general fiber of the momentum 
map of the KM systems is an affine part of a hyperelliptic Prym variety (or two 
copies of it), which can also be described as a hyperelliptic Jacobian. In order to 
describe precisely which affine part we determine the divisor which needs to be 
adjoined to each affine part in order to complete it into an Abelian variety. Since 
it is difficult to do this by using the maps $m we do this by performing Painleve 
analysis of the KM systems. 

The method that we use is based on the bijective correspondence between the 
principal balances of an integrable vector field (Laurent solutions depending on the 
maximal number of free parameters) and the irreducible components of the divisor 
which is missing from the fibers of the momentum map (see ||l[). 

We look for all Laurent solutions 

(20) a,(t) = lf;ap)i^ 



to the vector field (|13|) of the n-body KM system. The following lemma shows that 
any such Laurent solution of ( [l3| ) can have at most simple poles. We may suppose 
that r in ( pO| ) is maximal, i.e., af^ ^ for at least one i, and we call r the order 
of the Laurent solution. The order of pole (or zero) of ai{t) is denoted by r^, so 
r — maxi . 

Lemma 5.1. Let the Laurent series ai{t), i ~ 1, . . . ,n, given by be a solution 
to the vector field (j7^ of the n-body KM system. If at least one of the ai has a pole 
(for t — 0) then it is a Laurent solution of order 1. Moreover the orders of the pole 
( or zero ) of each ai {t) satisfy 

(21) n = 0^,-0!^,. 
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Proof. For s e N we find from (pO|): 

Res = I ^ ^ ' 

t=o a,(t) 1 0, s > 0. 

On the other hand, if we use ([l3| ) then we find 
We conclude that 



(22) 



(fc) _ (k) 



-Ti, /c = r — 1 
0, 0<k<r-2. 



Now substituting (pO|) into (13) and comparing the coefficient of the 



following equation [the indicial equation) is obtained: 

(23) -raf^^af )(«("), -al^,), z = l,...,n. 

If ai has a pole of order r > then a^^-* 7^ and ( p3| ) implies a^°\ — a^°\ = — r. 
Comparing with ( p2[ ) we see that we must have r = 1 and that ( pT| ) holds. 

□ 

Notice that in view of the periodicity of the indices (ai+„ = a^) the linear system 
l = (al+i~aS), i = 

has no solutions, so that at least one of the aj*^' vanishes. If, say, Oq*^' = 0^*^-^ = 
while af''^ ^ for i — 1, . . . , A: for some k in the range 1, . . . , 71 — 1 (this includes 
the case of a single i for which af^'' — 0) then the indicial equation specializes to 

«f = 1, 

a!+\-a|°\ = l, i = 2,...,k-l, 

which has no solution for k odd, and which has a unique solution (a^*^-*, . . . , a^°^) = 
(— /, 1, 1 — /, 2, . . . , —1, 1) for even k, k = 21. The other variables 0^*^-^ . . . , ai"-' can 
either be all zero, or they can constitute one or several other solutions of this type 
(with varying k = 21), separated by zeroes. Using periodicity the other solutions 
to the indicial equation are obtained by cyclic permutation. 

Thus we are led to the following combinatorial description of the solutions to 
the indicial equation of the n-body KM system. For a subset A of Z/n, and for 
p € Z/n let us denote by A{p) c Z/p the largest subset of A that contains p and 
that consists of consecutive elements (with the understanding that A{p) = when 
p ^ A). If we define 

I]„ ^ {Ac Z/n \ p e A^ #A{p) is even}, 

then we see that the solutions to the indicial equation are in one to one corre- 
spondence with the elements of I]„. In the sequel we freely use this bijection. For 
A S En we call the integer #A/2 its order, denoted by ord^. 
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For each solution to the indicial equation (i.e., for each A e I]„) we compute the 
eigenvalues of the Kowalevski matrix M, whose entries are given by 

dF,, 



do 



-{a 



(0) 



where Fi = ai{ai-i—ai+i), the i-th component of ([l3|). The number of non-negative 
integer eigenvalues of this matrix are precisely the number of free parameters of the 
family of Laurent solutions whose leading term is given by (aj"', . . . , ai"'') (see 0), 
hence we can deduce from it which strata of the Abelian variety, whose affine part 
appears as a fiber of the momentum map, are parameterized by it. 

Proposition 5.2. For a solution of the indicial equation corresponding to A E S„ 
the Kowalevski matrix M has n ~ ord A non-negative integer eigenvalues. 



Proof. In view of (21) the entries of M can be written in the form 

,(0) 



{l~ri)S^J, 



if a 







ifaf^O. 



Note also that, by using the Z/n action, we can assume that 1 G A, n ^ A, 
and that A is a disjoint union of A{pi), . . . , A{ps), with pi < p2 < ■ ■ ■ < Ps. Let 
k = Old A{pi). Then M has the following form 

/ 



M = 





El 




















E2 








D2 



I. \ 












Cs 


Es 


V 









I 

On the upper right corner the matrix has entry — Zi, and the blocks Ci,Di and Ei, 
i = 1, . . . , s, are matrices as follows: 

• Cj is a tridiagonal matrix of size 2li of the form: 

f I, \ 
1 -1 

1-h k-1 
2 



Q. = 



-2 



k-l 





-1 



l-k 


k 



1 



\ li J 

• Di is a diagonal matrix of the form Di ~ diag (1 + Z^, 1, . . . , 1, 1 + h+i), with 
the convention that if 1?^ is 1 x 1 then its only entry is 1 + li + Ij ; 



22 



RUI LOJA FERNANDES AND POL VANHAECKE 



• Ei is SL matrix with only one non-zero entry —li in the lower left corner; 
It is clear that the set of eigenvalues of M is the union of the set of eigenvalues 
of the Ci's and Di's. Now we have: 

Lemma 5.3. The eigenvalues of the matrix Ci are {±1, ±2, . . . , 

Assuming the lemma to hold we find that the number of negative eigenvalues of 
M is equal to X]i=i = Si=i ordA{pi) = ord^, so the proposition follows. 

So we are left with the proof of the lemma. We write I for h and we denote by ej 
the j-th vector of the standard basis of C^'. In the basis ei, 63, . . . , 621-3, 62/^1, 
e2Z, e2;_2, . . ■ , 64, 62 the matrix C; takes form 

A 
A 

where A is the transpose of the matrix 

/O 



A 







1 \ 








3 

i-i 



2 -1 
-2 



/ 



We show that this matrix has eigenvalues 1, —2, 3,.. .,(—!)' ^l. Then the result 
follows because the eigenvalues of C are ± the eigenvalues of A. 

For j I,..., I, let fj = [P-I, 2^-1, . . . , p-^]'^ and let Vj denote the span of 
fi, . . . , fj. For V — [vi, . . . , vi]'^ e C' we have that w S V, if and only if there exists 
a polynomial P of degree less than j such that Vk — P{k) ior k — 1, . . . , I. Since 
the fc-th component of Af,- is given by 

k{l-k+ ly-^ + (1 - k){l -k + 2)^-1 = (-l)J-ijF-i (^i + o(^^ 

we have that Af, C Vj , more precisely 

Af.ei-iy-'jt. + V,.,. 

This means that in terms of the basis {f,} the matrix A is upper triangular, with 
the integers 1, —2, 3, . . . , (— on the diagonal. 

□ 



' 1 free pa- 

-,(0)) ^ 



By the proposition above we can have a Laurent solution depending on n- 
rameters (a principal balance) only for the n choices of A given by {af\ . . . , 
(—1, 1, 0, . . . , 0) and their cyclic permutations. Let us check that these lead indeed 
to asymptotic expansions which formally solve (13). By §2 in Q|, these solutions 
are actually convergent and so they define convergent Laurent solutions. 



It suffices to do this for the solution {a\ 



(0) 



,(0)^ 



-1,1,0,. 



, 0) of the 



indicial equation. By (21) we know that the order of the singularities of this solution 
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(1,1,— 1,0,...,— 1) SO we have the foUowing ansatz for the formal 

ai{t) = - j+ai+ f3it + 0{t^), 
a2{t) = ^+a2 + P2t + 0{t^), 

aj{t) = + (3jt + 0[t^), 4<j<n-l. 
an(t)=l3nt + 0{e), 

If we replace these expansions in the equations (^) defining the n-body KM system 
we obtain the consistency equations: 

«! - "2 = 0, 

2/3l - /52 = -Q!lQ!2 - /?„, 

/3i - 2/32 = -aia2 + /Js, 

[ij = ajiaj^i - ttj+i), 4 < j < n - 1. 

They give exactly the n—1 free parameters ai, a^, . . . , a„_i, /Sa, The coefficients 
a*^*^^ ia[''\ . . . ,ali'^) for fc > 2 are then completely determined since they satisfy 
an equation of the form 

(M — kl) • a*-''-' — some polynomial in the ap' with j < k, 

and the eigenvalues of the Kowalevski matrix M are —1, 1, 2, by the proof above. 
This leads to the following result. 

Theorem 5.4. When n is odd the general fiber of the momentum map of the n-body 
KM system is an affine part of a hyperelliptic Prym variety, obtained by removing 
n translates of its theta divisor. When n is even the general fiber consists of two 
isomorphic components which admit the same description as in the odd case. In 
both cases the Prym variety admits an alternative description as a hyperelliptic 
Jacobian. 

6. Example: n=5 

In this section we study the 5-body KM system in more detail. Its phase space is 
four-dimensional and is given by /C5 = {(ai, 02, 03, 04, 05) | 0102030405 = 1}, with 
Lax operator 



/ 





Ol 








h-^ 


\ 




1 





02 















1 





03 















1 





04 




V 


hac, 








1 





/ 



The spectral curve det(a;Id —L) = is explicitly given by 

h+ \ = - Kx^ + Lx, 
h 



are (ri, . . . ,r„) = 
expansions: 



24 



RUI LOJA FERNANDES AND POL VANHAECKE 



where 

K ^ ai + a2 + as + ai + 05, 
L = aias + a2a4 + 0305 + a^ai + 0502- 

These functions are in involution with respect to the quadratic Poisson structure, 
given by {0^,0^} — (<5i.j+i — Si+ij)aiaj. It follows from the previous section that 
for generic k,l the afhne surface Vki defined by = A;, L = Z is an affine part of 
the Jacobian of the genus two Riemann surface Tr minus five translates of its theta 
divisor, which is isomorphic to r,-. As we have seen, an equation for Fr''' is given 

by 

(24) : y2 ^ (^3 - ku"^ + luf - Au. 

The two commuting Hamiltonian vector fields Xk and are given by 



hi 


= ai(a5 


- 02) 


O'l 


= 01(0305 


— O2O4) 


0.2 


= 02(01 


- 03) 


a'2 


= 02(0401 


- 03O5) 


03 


= 03(02 


- 04) 


03 


= 03(0502 


— O4O1) 


0,4 


= 04 (a3 


- 05) 


04 


= 04(0103 


- O5O2) 


0.5 


= 05(a4 


- Oi) 




= 05(0204 


- O1O3) 



The principal balance of Xk for which oi and 02 have a pole corresponds, ac- 
cording to Section |^, to the following solution of the indicial equations 

(or, of, of, of, of ) = (-!, 1,0, 0,0) 
and its first few terms are given by 

01 = -i + a - i(a2 + 2/3 + -f)t + 0{t^), 
a2 = j+a + i(a2 _ ^ _ 2^)t + 0{t^), 

(25) a3 = -/t + 0{t^), 

04 = (5 + 0(t2), 

05 =/3i + 0(<2). 

Here a, (3, 7 and 5 are the free parameters. If we look for Laurent solutions that 
correspond to the divisor to be added to Vu we find by substituting the above 
Laurent solution in K = k, L = I, 0102030405 = 1, 

2a + (5 = fc, 
2a(5 + /3-7 = Z, 
7/3(5 = -1, 

which means that the Laurent solution depends on two parameters /3 and 5, bound 
by the relation 

(26) ^k-S)S + P+j^=l, 

which is an (afhne) equation for the theta divisor, i.e., for r,-; it is easy to see that 
this curve is birational to the curve (p4| ). The other four principal balances are 
obtained by cyclic permutation from (|25|). 



HYPERELLIPTIC PRYM VARIETIES AND INTEGRABLE SYSTEMS 



25 



Vki can be embedded explicitly in projective space by using the functions with 
a pole of order at most 3 along one of the translates of the theta divisor and no 
other poles. Since the theta divisor defines a principal polarization on its Jacobian, 
the vector space of such functions has dimension 3^ = 9, giving an embedding in 
P*. One checks by direct computation that the following functions zq, . . . , zs form 
a basis for the space of functions with a pole of order at most 3 along the divisor 
associated with the Laurent solution (|2^) (the first two functions are obvious choices 
from the expression (^5|), while the others can be obtained from them by taking 
the derivative along the two flows): 





= 1, 


Zl 


= 0102, 


Zl 


= 010204, 


Zl. 


= 0102(01 + 05), 


Zi 


= 010204(03 + 04 + O5), 


Zb 


— 010204(01 — O2), 


Z6 


= 0i0204((03 + 04)01 — (04 - 


Z7 


2 2 
= 01020405, 


z% 


= 0i0204((04 + 05)^ + 0304) 



The corresponding embedding of the Jacobian in P** is then given explicitly on the 
affine surface Vu by (oi, . . . , 05) 1-^ (zq : ■ • ■ : zs)- By substituting the five principal 
balances in this embedding and letting i — > we find an embedding of the five 
curves Fi, . . . , r5 (in that order) which constitute the divisor Jac(rT-) \ Vu'- 



(0 : : : 1 : : 25 : : : -^3) 

(/?(52 : _/3252 . . _^2^3 : /35 : /3J : . Q . ;l _ ^^3) 

, ^ (1 : : /3(5 : : /3(5(fc-5) : /3(52 : -/3,5(/3 + 52 _fc5) : : /325(fc-<5)) 

{f3H -.0 : (36 : -f3S : P6{k - 6) : -f3S^ : 1 + f3S^ {6 - k) : 

: -6:-P5^{(3-(5-kf)) 

{PS^ : -S : : 6(6 - k) : P5 : -(35 : -(36^ : -1 : 1) 

The points on the divisor that correspond to the above Laurent solutions are the 
ones for which /? and S are finite; notice that all these points in P* are different. In 
order to determine the coordinates of the other points and the incidence relations 
between these points and the curves we choose a local parameter around each 
of the three points needed to complete ( ^6| ) into a compact Riemann surface: 

(a) 6 =l/u, (3 ^l/u^{l + 0{t)); 

(b) S=l/u, p^u^{l + 0{t)); 

(c) p S ^ ~u^{l + 0{t)). 
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Substituting these in the equations of the five embedded curves we find the following 
5 points (each one is found 3 times because it belongs to three of the curves Fj) 

pi = (0 : : : 1 : : : : : 0), 

J32 = (0 : : : : : : : : 1), 

P3 = (1:0:0:0:0:0:1:0: -k), 

P4 = (1 : : : : : : -1 : : 0), 

P5 = (0:0:0:0:0:0:0:1: -1). 

With this labeling of the points pi we have that Fj contains the points Pi-i,pi and 
Pi+i- As a corollary we find a 63 configuration on the Jacobian, where the incidence 
pattern of the 5 Painleve divisors and the 5 points pi is as in the following picture 
(to make the picture exact one has to identify the two points labeled ps, as well as 
the two points labeled p4 in such a way that the curves r2 and r4 are tangent, as 
well as the curves Fa and F5). 



Fs Fi F4 




F2 Fb 



Figure 2 

Obviously the order 5 automorphism 

(01,02,03,04,05) 1-^ (02,03,04,05,^1) 

preserves the affine surfaces Vki and maps every curve Fj and every point Pi to its 
neighbor. Since this automorphism does not have any fixed points it is a translation 
on Jac(F^), and since its order is 5 it is a translation over 1/5 of a period. Notice also 
that with the above labeling of points and divisors the intersection point between Fj 
and Fi+2 is Pi+i (so they are tangent), while the intersection points between Fj and 
Fj+i are Pi and pj+i . Dually, the divisors that pass through pi are precisely Fj_i , Fj 
and Fj+i. The usual Olympic rings are nothing but an asymmetric projection of 
this most beautiful Platonic configuration! 
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